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ABSTRACT 


As  the  target  radar  cross  section  (RCS)  continuously  decreases,  the  need  for  high- 
resolution  high-gain  radar  increases.  One  approach  to  high  resolution  is  to  use  distributed 
subarray  antennas  (DSAs),  because  of  limited  surface  available  on  many  radar  platforms. 

The  concept  of  distributed  subarray  antennas  is  examined  for  both  Multifunction 
Array  Radar  (MFAR)  and  Very  High  Frequency  (VHF)  applications.  By  combining 
distributed  subarrays  located  on  the  available  areas  of  a  constrained  platform,  the  MFAR 
and  VHF  DSA  can  achieve  increased  resolution  and  potential  reductions  in  cost  and 
complexity  compared  to  a  conventional  array.  The  two-way  pattern  design  of  DSA 
effectively  suppresses  the  undesired  grating  lobes  by  using  separate  transmit  and  receive 
antennas.  By  the  pattern  multiplication  principle,  the  grating  lobes  in  the  receive  pattern 
have  been  suppressed  by  proper  null  placement  of  subarray  in  the  receive  and  transmit 
patterns. 
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I.  INTRODUCTION 


A.  REQUIREMENTS  AND  OBJECTIVES 

The  platform  design  of  future  surface  combatants  has  changed  dramatically  with 
the  advent  of  stealth  technology  to  reduce  platform  signatures.  The  change  in  design 
philosophy  is  evident  in  the  proposed  new  DD-21,  as  illustrated  in  Figure  1.  The 
traditional  small-integrated  superstructure  of  vessels  allowed  few  areas  for  sensors  and 
weapon  systems.  On  the  other  hand,  the  threat  from  the  air  or  surface  is  ever  increasing 
such  that  the  performance  of  shipboard  radar  needs  to  provide  long-range  detection  and 
accurate  tracking.  This  implies  high  gain  and  physically  large  radar  antennas. 

Since  it  is  difficult  to  find  a  sufficient  area  for  a  large  array  onboard  a  ship,  it 
might  be  possible  to  use  several  smaller  noncontiguous  (separated)  areas  (subarrays)  and 
then  process  the  received  signal  coherently.  The  subarrays  may  be  far  apart  in  terms  of 
wavelength  and  therefore  grating  lobes  occur.  Grating  lobes  are  undesirable  because  they 
reduce  the  antenna  efficiency,  cause  ambiguities  in  angle  measurements  and  make  the 
radar  more  susceptible  to  jammers.  They  also  complicate  clutter  processing. 

The  main  objective  of  this  thesis  is  to  investigate  a  conceptual  design  of 
integrated  antenna  apertures,  which  are  composed  of  several  distributed  noncontiguous 
subarrays.  Together  the  entire  set  fonns  a  Multi-function  Array  Radar  (MFAR)  antenna 
that  provides  multi-tasking  with  high-resolution  high-speed  data  collection 
simultaneously.  This  thesis  will  show  the  trade-off  in  performance  and  cost  for  several 
possible  antenna  concepts. 

One  possible  approach  to  achieving  a  large  distributed  array  that  is  free  of  grating 
lobes  is  to  design  separate  receive  and  transmit  antennas.  For  example,  the  transmit 
antenna  might  consist  of  distributed  subarrays.  Then  a  separate  receive  antenna  would  be 
used  whose  pattern  has  nulls  in  the  direction  of  the  transmit  array’s  grating  lobes. 
Therefore  the  round  trip  or  “two  way”  gain  pattern  would  be  free  of  grating  lobes. 
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Figure  1.  The  Artist’s  concept  of  the  DD-21  [Ref.  1] 


This  approach  is  difficult  enough  for  narrow  band  radar.  However,  new  antenna 
designs  must  provide  wideband  perfonnance.  Furthermore,  due  to  the  limited  space,  it  is 
desired  that  the  antenna  serve  as  many  systems  as  possible:  several  radars, 
communications  systems,  and  electronic  warfare  systems.  Ideally,  the  antenna  should 
have  the  following  characteristics: 

1 .  Capability  to  execute  different  tasks  in  rapid  sequence 

All  sensors  and  electronic  devices  need  to  be  integrated  into  the  small  topside  of 
the  platform,  and  the  integrated  MFAR  should  have  surveillance,  tracking,  identification, 
lire  control,  missile  guidance  and  communication  capabilities.  The  bandwidth  of  the 
MFAR,  therefore,  must  be  wide  enough  to  provide  good  range  resolution  and  to  satisfy 
the  requirements  of  different  systems.  The  aperture  size  of  the  MFAR  must  be  electrically 
large  enough  to  give  line  angle  resolution  and  high  directivity.  The  side  lobes  of  the 
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antenna  pattern,  of  course,  have  to  be  low  enough  to  help  distinguish  the  target  from  the 
natural  or  man-made  clutter. 

2.  Capability  for  adaptive  digital  beamforming  (ADBF) 

Digital  beamforming  (which  requires  a  modulator  or  demodulator  at  each 
element)  is  preferred.  However,  the  large  number  of  array  elements  required  for  high 
resolution  reduces  the  possibilities  of  element-level  digital  beamforming  considerably. 
The  subarray-level  ADBF  could  provide  multiple  beams  and  electronic  scanning  to  carry 
out  the  multi-tasks  simultaneously  using  a  minimum  number  of  control  elements.  The 
MFAR  tracks  targets  accurately  with  low  sidelobe  sum  and  difference  beams,  which  are 
required  for  monopulse  tracking.  There  are  also  some  other  important  features  that  could 
be  provided  by  ADBF: 

(1)  Self-calibration  and  error  correction 

(2)  Adaptive  nulling  of  unwanted  interference  or  jamming  signals 

(3)  Spatial  control  of  the  radiated  energy  for  Low  Probability  of  Intercept 
(LPI)  operation 

3.  Affordability  with  acceptable  performance 

The  purpose  of  Multifunction  Subarray  Radar  (MFSAR)  is  to  fulfill  a  variety  of 
different  tasks.  Generally  the  cost  is  proportional  to  the  number  of  radiating  and  control 
elements,  as  well  the  complexity  of  construction  and  maintenance.  For  a  modern  high 
performance  radar  antenna,  the  basic  acceptable  antenna  parameters  are  the  following 
[Ref.  2]: 

Antenna  gains:  30  dB  (minimum) 

Target  tracking  accuracy:  5  mrad 
Azimuth  beam  steering:  ±45°  (Azimuth) 

Azimuth/elevation  coverage:  3607  60° 

Peak  sidelobe  level:  -45  dB 
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B.  PREVIOUS  APPROACHES 

There  are  two  basic  approaches  that  have  been  considered  as  a  solution  for  a 
wideband  multifunction  system  from  the  standpoint  of  the  array  aperture.  The  first  is  a 
segmented  approach,  where  each  segment  or  component  is  dedicated  to  a  separate  system 
and  bandwidth.  A  segmented  aperture  provides  better  performance  for  each  individual 
function  by  separating  the  size  and  position  of  the  antenna.  The  second  approach,  a 
shared  aperture,  however,  provides  multiple  frequency  operation  in  the  same  aperture. 
Generally  the  performance  at  a  single  frequency  is  not  optimum. 

The  proliferation  of  advanced  sensor  and  communication  systems  onboard 
military  platforms  has  led  to  a  multitude  of  systems.  A  U.S  Navy  Aegis  cruiser,  for 
instance,  has  over  100  antennas  and  the  number  is  expected  to  rise  as  new  systems  are 
added  [Ref.  3].  A  reduction  of  the  number  of  antennas  is  possible  using  the  shared 
aperture  concept. 

Ideally  a  shared  aperture  of  the  type  used  on  ships  should  operate  over  a  wide 
bandwidth,  for  example,  10  MHz  to  10  GHz  (three  decades).  If  the  aperture  is  shared  by 
separate  narrow  band  systems,  then  the  antenna  does  not  have  to  operate  over  the 
complete  continuous  range  of  frequencies,  but  only  at  the  operating  “sub-bands”  of  the 
component  systems.  These  types  of  antennas  are  referred  to  as  multi-band.  An  example  is 
the  Multifunction  Electronically  Scanned  Adaptive  Radar  (MESAR)  [Ref.  4]  that 
employs  a  dual  frequency  antenna.  The  work  in  MESAR  MFAR  programs,  which  began 
in  1977,  showed  significant  advantages  for  a  radar  operating  at  dual  optimal  frequencies. 
This  study  led  to  the  selection  of  1  and  10  GHZ  for  operating  frequencies  used  for 
surveillance  and  tracking  functions.  The  aperture  is  comprised  of  two  sets  of  antenna 
elements,  an  open-ended  waveguide  and  a  dipole,  for  the  low  and  high  frequency  bands, 
respectively.  The  proposed  system  will  therefore  perfonn  these  important  radar  functions 
at  maximum  efficiency. 
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C.  MULTIFUNCTION  SUBARRAY  RADAR 

Modern  phased  array  radars  are  used  in  a  multifunction  fashion,  with  the  required 
functions  being  search,  tracking  and  fire  control.  A  major  impetus  of  using  MFSAR  is 
cost  reduction.  For  example,  if  phase  control  is  at  the  subarray  level  and  each  subarray  is 
formed  by  5x5  elements,  then  the  total  number  of  transmit/receive  modules  is  reduced  by 
a  factor  of  25  for  contiguous  subarrays.  With  widely  spaced  subarrays,  not  is  there  a 
savings  in  control  elements,  but  also  higher  resolution.  When  the  subarrays  are  widely 
distributed  the  potential  resolution  increases  significantly.  The  disadvantage  is  the 
occurrence  of  grating  lobes  if  the  subarrays  are  widely  spaced  compared  to  the 
wavelength.  The  features  of  MFSAR  are  summarized  in  the  following  sections. 

1.  Advantages  of  MFSAR 

The  major  advantages  of  using  multifunction  subarray  radar  are  the  following: 

•  More  efficient  scheduling  of  the  track  and  search  functions  compared  to 
separate  systems  (track-while-scan,  TWS  capability). 

•  Rapid  steering  of  the  beam  to  the  desired  direction  as  needed. 

•  Formation  of  search  and  track  beams  with  maximum  flexibility. 

•  Savings  in  the  number  of  control  elements  and  decrease  in  the  complexity 
of  digital  beamforming. 

•  Increase  in  the  angular  resolution  by  distributing  subarrays  along  the 
superstructure  of  vessel. 

2.  Disadvantages  of  MFSAR 

The  major  disadvantages  of  using  multifunction  subarray  radar  are  the  following: 

•  Compromise  in  the  perfonnance  of  individual  functions  to  obtain  an 
optimum  balance  in  the  component  tasks. 

•  Extra  grating  lobes  introduced  by  large  subarray  spacing. 

•  Low  beam  efficiency  from  the  presence  of  grating  lobes. 
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D.  SCOPE  AND  ORGANIZATION 

1.  Scope 

In  this  research,  some  methods  of  correcting  the  above  disadvantages  are 
examined.  The  grating  lobe  problem  is  severe  for  subarray  spacing  over  a  wavelength. 
However,  it  is  possible  to  suppress  the  grating  lobes  through  control  of  the  subarray 
factor.  In  other  words,  the  nulls  of  the  subarray  pattern  are  placed  at  angles  where  the 
grating  lobes  occur.  There  is  a  limit  to  the  effectiveness  of  this  method,  as  the  gaps 
between  subarrays  get  larger.  Some  additional  improvement  is  achieved  if  separate 
transmit  and  receive  antennas  are  used,  and  grating  lobes  are  allowed  for  only  one  of  the 
two  antennas.  For  example,  if  the  transmit  antenna  has  grating  lobes,  then  the  receive 
antenna  will  not.  Furthermore,  the  nulls  of  the  receive  antenna  can  be  placed  at  transmit 
antenna  grating  lobe  locations.  Consequently  the  two-way  gain  pattern  will  have  no 
grating  lobes. 

An  additional  improvement  may  be  possible  with  digital  arrays  that  have 
amplitude  control  at  each  element.  The  transmit  array  factor  and  subarray  factor,  as  well 
as  the  receive  array,  can  be  phase  and  amplitude  weighted  to  give  low  sidelobe 
performance. 

This  thesis  demonstrates  that  distributed  antennas  of  MFSAR  can  provide  better 
performance  with  fewer  elements  than  the  conventional  MFAR  by  using  separated 
subarrays. 

2.  Primary  Research  Question 

There  are  two  related  research  questions  addressed  in  this  research: 

a)  how  to  increase  the  angular  resolution  of  shipboard  radar  through  wideband 
distributed  subarray  antennas,  and 

b)  how  to  suppress  the  grating  lobes  in  the  pattern  of  widely  spaced  subarrays  by  design 
of  the  subarray  pattern  or  design  of  a  two-way  pattern. 
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3.  Organization 

Chapter  II  provides  an  overview  of  the  array  theory  that  is  used  in  this  thesis.  A 
definition  of  the  array  factor  (AF)  is  given.  Antenna  parameters  such  as  beamwidth, 
directivity,  aperture  efficiency,  grating  lobes  and  pattern  scanning,  are  presented.  The 
assumptions  and  limitations  are  introduced. 

After  a  description  of  various  subarray  configurations,  Chapter  III  discusses  the 
different  methods  used  in  the  conceptual  design  of  an  MFSAR  antenna  and  the  concept  of 
digital  beamfonning  will  be  discussed  in  more  detail.  Some  ship  design  examples  are  also 
presented. 

Chapter  IV  provides  a  summary  and  conclusions,  followed  by  suggestions  for 
further  research  into  distributed  subarray  design. 

Appendix  A  provides  a  glossary  of  tenns  and  abbreviations  used  throughout  this 
thesis,  and  Appendix  B  provides  a  listing  of  MATLAB  codes  used  in  the  pattern 
calculations. 
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II.  SUMMARY  OF  ARRAY  THEORY 


This  chapter  discusses  the  fundamental  theory  of  phased  arrays.  The  material  can 
be  found  in  most  books  on  antennas,  such  as  [Ref.  5]  and  [Ref.  6]. 

A.  ARRAY  FACTOR  (AF) 

An  array  is  a  collection  of  smaller,  usually  identical  antenna  elements  that  are 
excited  with  complex  voltages  or  currents  to  obtain  a  desired  radiation  pattern.  For  most 
applications  the  elements  are  arranged  in  a  periodic  grid;  for  example,  in  two  dimensions 
it  could  be  a  rectangular  lattice.  The  array  can  be  divided  into  smaller  groups  of  elements 
called  subarrays  as  shown  in  Figure  2.  For  simplicity,  only  a  linear  array  is  shown,  and 
the  array  elements  are  assumed  isotropic. 


0=  90* 


*1 

< 


I 

>■< 


Isotropic  element 


2  =  (f 


Figure  2.  Example  of  an  array  of  contiguous  subarrays,  M  =  5  ,  N  =  3. 


9 


The  radiation  pattern  of  an  antenna  array  is  the  vector  sum  of  the  electric  field 
intensity  all  antenna  elements.  The  array  factor  for  N  uniform  subarrays,  each  composed 
of  M  isotropic  elements  along  the  x  -axis  is  given  by 


JV— 1 M- 1 


AF(0,</.)  =  YYl,Jme 


2n  l7t 

jn^-Dx(u-us)  jm — ax{u—us) 


n—0  m= 0 


M- 1  •  2/r 

v-1  Jm~r^x(u  us ) 

>  J  e 

/  j  m 

v-l  jn^j-O v  («-«s  ) 

2 

II 

o 

n—0 

(1) 


=  AFS  (&,</>)  x  AFc{9,  </>), 

where 

N,M  =  number  of  subarrays,  number  of  elements  in  a  subarray 
In,Jm  =  subarray,  element  weights  (excitation  currents  or  voltages) 
Dx,dx  =  subarray  center-to-center  spacing,  element  spacing 
u  =  sin  9  cos  (j) ,  direction  cosine  in  spherical  coordinates 
us  =  sin  9S  cos  (f)s  ( 9S ,  (j)s  is  the  scan  direction) 


AFs  =  subarray  pattern  detennined  by  the  arrangement  of  elements  in  a  subarray 
AFc  =  configuration  pattern  determined  by  the  arrangement  of  subarrays. 

Equation  (1)  assumes  that  all  subarrays  are  identical  and  the  elements  have 
uniform  spacing.  The  pattern  of  the  array  is  then  the  product  of  the  element  factor  and  the 
array  factor.  Isotropic  elements  have  been  assumed  and,  therefore,  the  element  factor 
(EF)  is  1 .  The  use  of  subarrays  allows  the  array  factor  to  be  separated  based  on  the 
geometry  and  excitation  of  subarrays. 
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B.  ANTENNA  PARAMETERS 

1.  Beamwidth  Between  First  Nulls  (BWFN) 

The  beamwidth  between  the  first  nulls  is  defined  by  the  subtended  angle  of  the 
mainlobe  in  the  antenna  pattern.  If  the  element  and  subarray  weights  in  Equation  (1)  are 
assumed  unity,  then  the  normalized  antenna  pattern  is  given  by: 


where 


\AFnorMt  = 


sin(M)  sin(^) 


M sin(^)  N sin(^) 


£ 


l7T 

T 


dx(u  —  us) 


(2) 


y/=—Dx(u-us). 

The  subscript  “norm”  denotes  the  normalized  array  factor.  The  BWFN  is  determined  by 
the  first  zero  of  the  numerator.  Since  the  phase  tenn  of  the  subarray  is  changing  much 
faster  than  the  phase  of  elements  ( NDx  »  Mdx ),  the  BWFN  is  detennined  mainly  by  the 

phase  of  subarrays  given  by  u-us  =  ±1  / NDX  so  0BWFN  =  2sin_1(l  / NDx)  in  the  plane 
(/)  =  0  .  The  BWFN,  therefore,  is  determined  by  the  entire  aperture  size  (i.e.,  the  distance 
between  the  edges  of  the  two  farthest  separated  subarrays,  L  ~  ND  ). 

2.  Half  Power  Beamwidth  (HPBW) 

The  half-power  beamwidth  of  a  linear  array  is  defined  as  the  angular  separation 
between  two  points,  one  on  each  side  of  the  main  beam  maximum,  at  which  the  power  is 
reduced  by  half.  The  HPBW  is  determined  mainly  by  the  subarray  configuration  pattern. 
For  a  half-wavelength  spaced  linear  array  with  unifonn  excitation,  the  HPBW  can  be 
approximated  by  0HPBW  ~  50.8/1/ L  in  degrees,  where  L  ~  NDx  is  the  entire  length  of  the 
aperture. 
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3.  Directivity  and  Gain 

Directivity  is  defined  as  the  ratio  of  the  maximum  radiation  intensity  in  the  main 
beam  to  the  average  radiation  intensity  [Ref.  7].  For  a  two-dimensional  unifonn 
rectangular  array  of  isotropic  elements,  the  directivity  can  be  approximated  by 
D  =  4 7TT]LxLy  cos#  /  A2 ,  where  ij  is  the  aperture  efficiency,  Lx,Ly  are  the  dimensions  of 

the  array  in  th ex,  v  directions,  and  0,  is  the  scan  direction.  This  approximation  is  good  as 
long  as  there  are  no  grating  lobes  in  the  visible  region.  Subsequently  it  is  assumed  that 
there  are  no  other  losses  other  than  those  due  to  amplitude  tapering.  In  this  case  the 
directivity  and  gain  are  equal. 

4.  Aperture  Efficiency 

Aperture  efficiency  is  a  measure  of  how  efficiently  the  antenna  physical  area  is 
utilized.  If  the  element  or  subarray  amplitudes  are  not  equal,  the  aperture  efficiency  is 
given  by  the  taper  efficiency.  For  a  planar  array  in  the  x  —  y  plane 


?I  = 


Nx  Ny 

SIX 


Nx  Ny 


a'aEXIc 


n- 1  m=l 


(3) 


where  Nx  and  Ny  are  the  number  of  elements  in  the  x  and  y  directions,  respectively,  and 
Imn  is  the  amplitude  of  the  mn  th  element  weight.  The  weighting  can  be  applied  to  either 
the  AFc  or  AFs  summation. 

5.  Grating  Lobes 

Grating  lobes  occur  when  more  than  one  period  of  the  array  factor  appears  in  the 
visible  region  (±90°);  when  either  the  element  or  subarray  spacing  in  the  Equation  (1)  is 
more  than  one  wavelength.  For  MFAR  it  is  likely  the  subarray  spacing  is  larger  than  one 
wavelength,  and  thus  grating  lobes  will  exist,  as  shown  in  Figure  3.  In  theory,  the  grating 
lobes  could  be  suppressed  by  using  either  unequally  sized  or  randomly  spaced  subarrays, 
at  the  expense  of  higher  complexity,  average  sidelobe  level  increase,  and  difficulties  in 
applying  ADBF  [Ref.  8].  Another  possible  means  of  suppressing  grating  lobes  is  to  use 
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two  separate  antennas  for  transmit  and  receive  (i.e.,  a  quasi-monostatic  radar).  The  two- 
way  gain  GtGr  occurs  in  the  radar  range  equation  for  the  signal-to-noise  ratio  (SNR) 


SNR  = 


P.G.G^cj 
(4 7rfR4KTsBn  ’ 


(4) 


where  Pr  is  the  received  target  signal,  Pt  is  the  transmit  power,  a  the  target  RCS,  X  the 
wavelength  and  R  the  range.  The  product  KTsBn  is  the  noise  power,  K  is  Boltzmann 
constant,  Ts  system  noise  temperature  and  Bn  radar  bandwidth.  Grating  lobes  can  be 

allowed  in  the  receive  pattern  and  then  eliminated  by  placing  transmit  nulls  in  those 
directions.  The  cost  is  usually  a  reduction  in  transmit  antenna  efficiency. 


Figure  3.  Grating  lobes  produced  by  subarrays  Dx  =  5A  (dx  =  0.5/1  ,N  =  5,M  =  5) . 
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6.  Electronic  Scanning 

Electronic  scanning  is  defined  as  a  method  of  positioning  a  beam  with  the  antenna 
aperture  remaining  fixed.  The  basic  electronic  scanning  techniques  are  phase  shifting, 
true  time  delay,  frequency  scanning,  and  feed  switching.  Some  possible  control 
configurations  for  beamforming  are  shown  in  Figures  4  through  6.  To  avoid  the  problems 
of  beam  squinting  and  broadening  that  occur  over  wide  frequency  bands,  it  is  desirable  to 
have  true  time  delay  at  each  element.  This  is  difficult  to  achieve  using  conventional 
microwave  beamforming.  However,  it  has  been  accomplished  using  photonic  devices 
(optical  fiber)  and  through  signal  processing  in  digital  antennas.  Subsequently,  true  time 
delay  is  assumed  at  each  element. 


PHASE 

SHIFTERS 


SUBARRAYS 

IN/OUT 


TRUE  TIME 
DELAY 


Figure  4.  True  time  delay  at  the  subarray  level  with  phase  shifters  at  the  element 
level. 


14 


Figure  6.  Digital  beamforming. 
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c. 


RADAR  SYSTEM  CONSIDERATIONS 


The  radar  range  equation  (4)  gives  the  fundamental  relationship  between  radar 
design  parameters.  A  main  objective  of  this  thesis  is  to  propose  a  possible  conceptual 
design  of  distributed  subarray  antennas,  using  separate  transmitter  and  receiver  antennas, 
which  is  a  fonn  of  quasi-monostatic  radar. 

1.  Maximum  Detection  Range 

From  Equation  (4),  the  maximum  detection  range  of  the  radar  for  a  given  signal- 
to-noise  ratio  can  be  rewritten  as 

R  4  _  P,GtG,X(J  _  PfrAj; 
max  W  (SNR),,  KTsBn  {4nf  KTsBn  (SNR)n 

Pg,sGAer<y  (5) 

(4x)2 KTsFn(BT)fp(SNR)H  ’ 

where 

P  =  average  transmitter  power 
Aer  =  receiver  effective  aperture  area 


Fn  =  receiver  noise  figure 

B  =  receiver  bandwidth 
T  =  pulse  width 
/  =  pulse  repetition  frequency 

( SNR)n  =  minimum  signal-to-noise  ratio  after  n  -pulses  are  integrated. 

This  equation  illustrates  several  important  tradeoffs  in  the  design  of  radar 

(1)  Power- aperture  product:  One  of  the  most  important  measures  of  the 
capability  of  long-range  surveillance  radar  is  the  product  of  average  power 
and  effective  aperture  size  (PngAer) ,  which  controls  the  amount  of  power 

transmitted  by  the  aperture.  If  the  aperture  size  could  be  larger,  the  average 
power  required  for  detecting  the  target  at  a  given  range  would  be  decreased 
when  all  other  factors  are  constant. 
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(2)  Frequency  dependence:  Although  the  frequency  does  not  appear  explicitly,  it 
is  easier  to  achieve  high  power  at  low  frequencies  because  high  voltages  are 
applied  to  larger  distances  and  breakdown  can  be  avoided.  Furthermore, 
atmospheric  loss  is  less  at  lower  frequencies. 

2.  Angular  and  Range  Accuracy 

Range  and  angular  resolution  provide  improved  accuracy,  since  accuracy  is 
inversely  proportional  to  signal  bandwidth  and  directly  proportional  to  beamwidth  as 
[Ref.  9] 


<?e  = 


e, 


HPBW 


K,j2  SNR, 


(6) 


where 


<J£  =  angular  accuracy 


0HPBW  =  antenna  beamwidth 


Kg  =  gradient  of  the  difference  beam  in  the  monopulse  antenna  configuration. 

Once  again,  the  angular  accuracy  is  proportional  to  the  beamwidth  0HPBW  and  inversely 
proportional  to  square  root  of  SNR  in  this  equation. 

Since  the  range  resolution  is  A R  =  ct ! 2  =  c/2B  [Ref.  3],  the  range  accuracy  is 
improved  by  having  a  wide  instantaneous  bandwidth. 
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D.  SUMMARY 

This  chapter  has  presented  the  basic  theory  of  arrays  and  discussed  some 
fundamental  radar  system  design  tradeoffs.  This  research  is  primarily  concerned  with 
increasing  radar  resolution  by  an  increase  in  the  antenna  aperture.  Because  large 
unobstructed  smooth  surfaces  are  limited  on  ships  and  aircraft,  it  is  necessary  to  construct 
large  arrays  by  combining  signals  from  distributed  subarrays.  Unfortunately,  this  results 
in  grating  lobes. 

In  order  to  suppress  the  grating  lobes  the  subarray  sizes,  locations  and  excitations 
can  be  adjusted.  Furthermore,  separate  transmit  and  receive  antennas  can  be  employed, 
and  designed  so  that  the  two-way  pattern  has  suppressed  grating  lobes.  These  last  two 
approaches  are  examined  in  subsequent  chapters.  To  simplify  the  analysis  the  following 
assumptions  are  made: 

1 .  The  element  factor  is  neglected, 

2.  Generally,  a  linear  array  is  used  with  x  being  the  array  axis, 

3.  The  antenna  is  at  the  center  of  a  spherical  coordinate  system,  where  the 
x  —  y  plane  is  the  earth’s  surface  and  z  the  zenith  direction, 

4.  Mutual  coupling  is  neglected,  and 

5.  A  constant  frequency  or  time-harmonic  wave  is  assumed  ( eJwr  time 
dependence),  and  therefore  phaser  quantities  appear  in  the  equations. 
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III.  DISTRIBUTED  SUBARRAY  ANTENNA 


A.  SUBARRAY  METHODS  AND  CONFIGURATIONS 

The  concept  of  subarraying  arises  from  the  requirement  of  modem  radar  for  high 
resolution.  When  the  targets  to  be  tracked  have  very  small  RCS,  and  the  radar  must  have 
fine  angle  resolution  to  track  multiple  closely  spaced  targets,  the  antenna  apertures  have 
to  be  electrically  large.1  Because  of  the  large  number  of  array  elements  required, 
grouping  of  elements  to  reduce  the  beamfonning  complexity  and  control  cost  is 
inevitable.  The  subarrays  can  be  divided  into  several  types  as  shown  in  Figure  7.  They 
are  discussed  in  the  following  section. 

1.  Subarray  Configurations 

In  terms  of  control  there  are  generally  four  types  of  subarray  configurations,  as 
shown  in  Figure  7: 

1)  Amplitude  and  phase  control  at  each  element  (Figure  7(a))  —  this  allows 
complete  control  of  both  AFs  and  AFc , 

2)  Phase  only  control  at  each  element  (Figure  7(b))  —  this  allows  scanning  of 
both  AFc  and  AFs , 

3)  Amplitude  and  phase  control  of  the  subarray  pattern  AFc  (Figure  7(c))  and 

4)  Amplitude  and  phase  control  of  the  array  factor  AFc  (Figure  7(d)). 


1  Synthetic  aperture  radar  (SAR)  is  another  solution  to  the  resolution  problem,  but  it  has  its  own 
disadvantages  and  limitations  [Ref.  10]. 
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AMPLITUDE/PHASE  (a) 


PHASE-ONLY  (b) 


SELECT  ELEMENTS  (c)  SUB  ARRAY  WEIGHTING  (d) 

Figure  7.  Typical  subarray  control  configuration  [From  Ref.  1 1] 


Amplitude  and  phase  control  at  each  element  is  undoubtedly  the  most  costly  but 
also  the  most  desired  in  terms  of  beam  pointing  and  sidelobe  level  (SLL).  The  second 
alternative  is  to  control  element  phase  only,  which  is  very  attractive  since  in  a 
conventional  phased  array  the  required  controls  are  available  at  no  extra  cost.  However, 
the  grating  lobes  in  the  DSA  approach  severely  distort  the  pattern,  and  a  phase-only 
correction  is  not  very  effective.  In  this  thesis,  it  is  assumed  that  complex  weights  can  be 
applied  at  both  the  element  level  and  the  subarray  level.  Furthermore,  the  weights  have 
true  time  delay  behavior.2 


2  As  a  practical  matter,  if  complex  weights  can  be  applied  at  the  element  level,  then  there  is  no  need 
for  weighting  at  the  subarray  level. 
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2.  Basic  Properties  of  DSA  Pattern 

Due  to  the  large  spacing  between  subarrays  in  a  periodic  DSA,  there  are  many 
grating  lobes  in  the  visible  region.  For  example,  if  the  spacing  between  subarrays  is  5 
wavelengths,  there  will  be  10  grating  lobes  in  the  [-;r/2,;r/2]  visible  region,  that  are 

located  at  the  angles  sin_1(±/>/5)  radians  (p  =  1,2,3,  — )  as  shown  in  Figure  3.  The 
following  are  some  basic  characteristics  of  the  periodic  subarray  configuration: 

•  The  array  pattern  of  identical  periodic  subarrays  can  be  separated  as  a 
multiplication  of  a  subarray  factor  ( A  Fs )  and  the  subarray  geometry  and 

excitation  factor  (AFc)  as  shown  previously  in  Equation  (1), 

•  The  number  and  the  intensity  level  of  sidelobes  between  two  adjacent 
grating  lobes  are  dependent  on  the  number  and  position  of  the  subarrays. 
Similar  to  the  uniform  linear  array,  there  will  be  N  -  2  sidelobes  with  a 
peak  intensity  level  about -13  dB  in  the  unity  weighting  case,  and 

•  The  number  of  subarrays  does  not  have  any  effect  on  the  number  of  the 
grating  lobes  as  long  as  subarrays  are  spaced  approximately  equal.  The 
intensity  level  of  the  sidelobes  is  decreased  with  an  increase  in  the  number 
of  subarrays.  This  implies  that  having  multiple  small  subarrays  that  are 
spaced  uniformly  gives  better  performance  in  both  SLL  and  grating  lobes 
than  having  two  large  subarrays  spaced  far  apart. 

3.  Methods  for  DSA  Pattern  Design 

There  are  two  approaches  to  reducing  grating  lobes.  The  first  is  based  on  pattern 
multiplication,  as  illustrated  in  Equation  (1).  The  grating  lobes  in  one  factor  can  be 
suppressed  by  placing  nulls  of  the  other  factor  coincident  with  the  grating  lobes.  The  key 
to  this  approach  is  that  the  grating  lobes  and  nulls  are  periodic  in  arrangement.  A  second 
approach  to  reducing  grating  lobes  is  to  use  an  irregular  spacing  or  unequal  sized 
subarrays,  thereby  reducing  the  peak  grating  lobes  by  redistributing  the  energy  into  the 
sidelobe  regions. 
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The  process  of  pattern  design  can  be  divided  into  two  stages.  The  first  stage  is  to 
lower  the  SLL  and  eliminate  the  grating  lobes  as  much  as  possible  by  synthesis  of  the 
appropriate  AFc  and  AFs  functions.  For  the  subarray  spacings  of  interest  the  first  pair  of 

grating  lobes  (one  on  either  side  of  the  main  lobe)  are  not  easy  to  suppress  because  of 
their  narrow  width  and  close  proximity  to  the  main  lobe.  The  second  step  is  to  suppress 
the  remaining  grating  lobes  by  either  a  high  contrast  (i.e.,  low  sidelobes  relative  to  the 
mainlobe)  transmit  pattern  or  by  specific  placement  of  the  nulls  in  the  transmit  pattern. 

There  are  several  simple  methods  used  in  the  pattern  synthesis  of  DSA  in  this 
thesis.  The  process  of  direct  nulling  by  pattern  synthesis  is  discussed  next. 

a.  Pattern  Multiplication  Principle 

Since  the  antenna  is  composed  of  periodic  subarrays  with  the  same 
spacing,  their  pattern  could  be  predicted  from  Equation  (1)  and  the  null  locations 
manipulated.  The  angular  directions  of  grating  lobes  in  the  configuration  pattern  AFC  are 
given  by 


u  =  us  ± 


(7) 


where  p  is  a  non-zero  integer.  The  nulls  of  the  uniform  subarray  pattern  AFs  are  located 
at 


u  =  us  ± 


p/ 1 


(8) 


The  condition  for  suppressing  grating  lobes  is  DJdx=M.  Unfortunately  this  is  exactly 

the  condition  of  contiguous  subarray,  which  contradicts  our  DSA  approach.  Since  the 
nulls  of  the  uniform  subarray  pattern  are  located  in  2n  /  M  angular  increments,  if  the 
element  number  M  is  chosen  correctly,  some  ratio  of  grating  lobes  in  the  configuration 
pattern  will  be  suppressed.  For  example,  if  Dx  / dx  =  5A/0.5A ,  and  let  M  =  5 ,  the  even 
integer  numbers  of  grating  lobes  will  be  suppressed  as  shown  in  Figure  8  (a). 
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This  only  solves  half  of  the  problem.  The  two-way  pattern  synthesis  gives 
additional  freedom  to  use  the  nulls  of  the  transmit  pattern  to  suppress  the  remaining 
grating  lobes  in  the  receive  pattern.  The  two-way  pattern  is  defined  as 

AFlway  =  AFTx  x  AFRx  =  AFTx  x  AFs  x  AFc .  (9) 

For  example,  the  receive  pattern  of  16  weighted  subarrays,  and  transmit  pattern  of  20 
weighted  elements,  spaced  in  0.5/1  is  shown  in  Figure  9.  In  this  example  the  transmit 
antenna  has  a  Chebyshev  distribution. 


23 


Receive  pattern  of  16  subarrays  spacing  5X 


Pattern  of  single  subarray  which  is  composed  by  5  elements 


Pattern  of  16  subarray  configuration 


(b) 


(c) 


Figure  8.  The  grating  lobes  suppressed  by  5  element  subarrays  pattern,  (a)  total 
receive  array  pattern,  (b)  subarray  factor  AFs ,  and  (c)  subarray  configuration  factor  AFc . 
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Two  Way  pattern  of  1 6  subarrays 


(a) 


(b) 


(c) 


Figure  9.  Example  of  pattern  multiplication  principle  (a)  two-way  pattern  AF2  , 
(b)  receive  pattern  AFRx ,  (c)  transmit  pattern  AFTx . 
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b.  Modified  Minimax  Algorithm  to  Find  the  Geometry  of  Lowest 
Sidelobes  by  Perturbation  of  Subarray  Location 

There  are  many  different  approaches  to  the  design  of  low  sidelobe  patterns 
in  phased  array  systems.  Since  the  goal  of  the  first  step  is  to  minimize  the  sidelobes  and 
grating  lobes,  one  such  approach  is  to  displace  the  subarray  locations  based  on  the  Dolph- 
Chebyshev  theorem  [Ref.  12]. 

Every  subarray  considered  is  identical  (same  AFs )  and  has  symmetric 
weights  about  the  center  of  the  array.  The  AFC  can  be  represented  as 


where 


K~l  2k  D 

AFc(u )  =  ^ a k  cos( - -u) ,  0  <  u  =  sin#  <  1, 

X 


k= 0 


(10) 


N 

K  =  —  when  N  even, 
2 


jV  + 1 
2 


when  N  odd,  N  is  the  number  of  subarrays  and 


Dk  =  distance  between  the  Mi  subarray  and  0th  subarray. 

For  odd  length  arrays  the  weights  to  be  applied  on  the  individual  subarrays  are  found 
from  L  =  a ,  1=1=  a„  /  2  and  in  the  even  case  I  =  I,  =  al  2 . 

By  fixing  the  location  of  the  center  and  edge  subarrays,  the  lower  bound 
on  the  beamwidth  of  the  main  lobe  of  AFc  is  approximately  u0  ~  0.886A/  L  where  L  is 

the  length  of  aperture,  as  the  initial  guess.  The  Minimax  algorithm  is  used  by  setting  the 
lower  and  upper  bound  for  location  of  other  subarrays,  and  the  obtaining  subarray- 
optimized  positions  by  minimaxing  the  Chebyshev  basis  function  in  the  region  u0  <  u  <  1 . 
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The  application  of  this  method  to  the  16  uniform  subarrays  of  5  unifonn 
elements  each  is  shown  in  Figure  10.  Compared  to  the  original  receive  pattern  of  16 
Hamming-weighted  subarrays  with  equal  spacing  in  Figure  8(a),  although  the  peak 
grating  sidelobes  are  decreased  from  4  dB  to  11  dB,  the  average  sidelobe  level  becomes 
much  higher  because  of  the  perturbation  in  locations.  Also  note  that  the  narrow 
beamwidth  is  similar  to  the  uniform  weighted  subarray.  The  locations  of  the  subarrays  in 
wavelengths  are  Dn=  [0  2.5000  11.6208  15.1714  19.2073  24.6143  29.2136 

35.0000  40.0000  45.7864  50.3857  55.7927  59.8286  63.3792  72.5000 

75.0000]. 


Figure  10.  One-way  array  pattern  of  16  location-perturbed  unifonn  subarrays 
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The  same  Minimax  algorithm  can  also  be  extended  to  synthesis  of  the 
weights  of  the  desired  transmit  pattern.  It  will  find  the  suitable  weights  of  the  transmit 
elements  for  selective  nulling  angles  and  the  specified  sidelobe  level,  which  match  the 
grating  angles  in  the  receive  pattern.  The  objective  function  becomes  [Ref.  13] 

min { max H \W(uk ){D(uk)  ~P(uk))\}  =  8  (11) 

where 

u  =7T  sin#  ,  the  set  of  spatial  frequencies  the  response  is  optimized 
W  =  the  error  weighting 
D  =  the  desired  response  (intensity  level) 

8  =  deviation  between  the  pattern  and  desired  response 

P(uk )  =  Yjak  cos (&u) ,  4  =  — k- . 

k= 0 

This  theorem  could  be  formulated  as  a  matrix  equation  A  a  -  D  where 


COS 

cos  £ku0 

8 

W(u0) 

A  = 

cos  ^0u{ 

COS^Wj 

-8 

W(ux) 

cos  <^QuK+1 

COS  ^KUK+l 

(-1  )K+18 
m<M)_ 

ll 

1  $ 

8]t 

D(u0 )  D(ux)  ■ 

••  D(uk) 

The  coefficients  ak  that  form  the  solution  to  this  system  can  then  be  used  as  the  element 

weights  of  the  transmit  array.  This  approach  was  tried  but  the  null  depths  were  not  always 
sufficient  to  suppress  the  grating  lobes. 
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c.  Modified  Weighting  Method  for  Distributed  Subarray 

The  conventional  weighting  distributions  can  be  applied  to  nonuniformly 
spaced  arrays  by  sampling  the  continuous  distribution  at  the  appropriate  points  that 
correspond  to  the  element  locations.  The  phase-shift  of  elements  for  scanning  needs  to  be 
modified  accordingly  [Ref.  14]. 

Let  the  number  of  the  subarrays  be  N,  the  number  of  array  elements  in 
the  n  th  subarray  M  ( n )  =  M  (identical  subarrays),  the  distance  between  the  nth  and  Oth 
subarrays  are  L(n)  wavelengths,  and  0k  the  angle  from  nonnal  direction  of  the  array. 
Then  the  phase-shift  of  the  mi h  element  in  nth  subarray  is 

2  7r 

(Pmn  =L(n)—sxndk  +m—dxsmdk.  (12) 

The  modified  Hamming  weighting  function  for  a  nonunifonn  spaced  subarray  is  written 
as 


where 


WH(n)  =  0.54  -0.46  cos[^^], 

A/ 


(13) 


ld{n)  =  L(n)/ dx  nth  subarray  position  relative  to  the  element  spacing 


Ld  =  Lx  /  dx  aperture  dimension  relative  to  the  element  spacing. 

A  similar  result  could  be  achieved  if  this  modified  weighting  method  is  applied  to  the 
elements  in  all  subarrays,  which  yields  a  smoother  sampled  distribution  than  weighting  at 
the  subarray  level  only.  However,  the  array  factor  cannot  be  separated  because  of 
different  subarray  factors  AFS.  Also  ld(n)  has  to  be  changed  to  the  element  relative 

position  ld(mn)  in  the  aperture  to  apply  the  weighting  function  of  Equation  (13). 
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Figure  1 1  shows  the  receive  pattern  of  modified  Hamming  weight 
application  on  the  elements  of  the  same  subarrays.  Compared  to  Figure  9(b),  there  is 
more  suppression  of  the  grating  sidelobes  at  the  angles  #  =  ±23°, ±53°,  but  almost  the 
same  level  everywhere  else.  The  patterns  for  weighted  subarrays  versus  weighted 
elements  do  not  differ  too  much  in  beamwidth  and  sidelobe  level.  The  subarray  level 
control  scheme  is  more  desired  because  of  its  simplicity. 

For  reference,  Figure  12  shows  the  pattern  of  16  uniform  weighted 
subarrays.  The  pattern  of  the  Hamming  weighted  elements  or  subarrays  is  “cleaner”  in 
terms  of  sidelobe  level  especially  in  the  proximity  of  the  mainlobe.  The  Hamming  pattern 
has  a  slightly  wider  beamwidth. 


Figure  11.  The  DSA  receive  pattern  with  modified  Hamming  weights  on  array 
elements 
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Figure  12.  The  DSA  Receive  pattern  with  uniform  weighted  subarray  elements 
( N  =  1 6,  M  =  5,  Dx  =  5  A,  dx  =  0.5  A ) 

d.  Equiripple  (Parks-McClellan)  Design  Method  of  High  Contrast 
Transmit  Pattern 

Linear  antenna  arrays  are  in  many  ways  analogous  to  one-dimensional 
digital  filters.  Restriction  of  the  pattern  synthesis  problem  to  that  of  discrete  arrays  of 
finite  spatial  extent  makes  the  problem  similar  to  that  of  finite  impulse  response  (FIR) 
digital  filters.  When  the  ideal  time  delays  are  used,  wideband  pattern  synthesis  reduces  to 
the  narrow  band  case,  with  each  element’s  delayed  waveform  receiving  a  single  real 
weight.  In  classic  narrow  band  pattern  synthesis,  an  equiripple  weighting  with  narrowest 
beamwidth  for  a  given  sidelobe  level  was  proposed  in  [Ref.  15]. 
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The  Parks-McClellan  algorithm  is  based  on  an  iterative  algorithm,  which 
minimizes  the  maximum  amplitude  of  the  ripple  (side  lobes)  present.  By  the  specified 
angle  of  the  passband  (BWFN),  stopband  (spatial  directions  outside  the  mainlobe), 
frequency  response  (intensity  level)  and  the  maximum  deviation  from  the  frequency 
response  the  weighted  least  squares  algorithm  (WLS)  provides  an  optimal  approximation 
to  the  desired  pattern  in  the  least  squares  sense.  This  ensures  that  the  power  present  in  the 
stopband  will  be  a  minimum.  The  calculation  of  an  array  factor  of  35  linear  elements 
spaced  XU  using  MATLAB’s  built-in  Remez  function  is  shown  in  Figure  13,  and  the 
applied  weights  for  each  element  are  shown  on  Figure  14. 


Figure  13.  Transmit  pattern  of  35  equiripple-weighted  elements 
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Figure  14.  Weights  of  the  transmit  elements 

B.  POSSIBLE  SHIPBOARD  MFAR  DSA  DESIGN 

There  are  several  reasons  why  a  distributed  subarray  antenna  might  be  used  on  a 
ship.  One  is  the  limited  surface  area  available  for  antenna  placement.  At  low  frequencies, 
the  open  surface  areas  are  small  compared  to  the  wavelength,  and  narrow  beamwidths  are 
not  possible  with  a  single  contiguous  array.  By  adding  more  subarrays  on  other  areas  of 
the  ship,  they  can  be  processed  as  a  DSA  resulting  in  a  half  power  beamwidth  of 
approximately  XI L  radians.  Any  combination  of  transmit  and  receive  functions  could  be 
used.  In  Figure  15,  for  example,  1  and  2  transmit  only;  3  receives  only,  etc.  This 
approach  can  be  applied  to  frequencies  where  the  subarray  spacing  is  in  the  range  of  1  to 
5  wavelengths.  For  an  Aegis-sized  cruiser  this  would  be  in  the  VFIF  to  UHF  frequency 
regions. 

At  higher  frequencies  a  DSA  might  be  used  to  reduce  cost  and  weight.  As 
illustrated  in  Figure  16  the  original  aperture  size  is  large  enough  (in  wavelengths)  to 
provide  a  sufficiently  narrow  beam,  X!  L.  Some  weight  can  be  eliminated  by  removing 
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selected  areas  of  the  array,  as  shown  in  Figure  17.  The  thinned  areas  between  the 
subarrays  can  be  used  for  other  sensors;  for  example,  radar  or  communication  antennas  at 
other  frequencies.  The  penality  is  grating  lobes,  although  they  can  be  suppressed  using 
the  techniques  described  previously. 


34 


1.  MFAR  DSA  Configuration  by  Thinning 

One  possible  antenna  configuration  for  a  shipboard  DSA  is  shown  in  the  Figure 
17.  C-band  is  chosen  as  the  operating  frequency  for  this  design.  The  transmit  array  is 
composed  of  the  center  35x35  uniform  elements  (yellow),  which  also  can  be  used  as 
part  of  the  receive  antenna  (blue  squares).  In  theory  the  transmit  pattern  could  use  any 
number  of  elements,  with  the  performance  improving  with  the  number  of  elements 
(narrower  transmit  beamwidth  and  higher  gain). 

The  receive  antenna  uses  digital  beamforming  with  16x8  subarrays,  each 
comprised  of  5x5  isotropic  elements.  Each  subarray  can  have  any  number  of 
independent  functions  such  as  communication,  missile  guidance,  sidelobe  cancellation, 
etc.,  and  they  are  grouped  independently  for  each  function.  The  following  section  shows 
patterns  for  16  by  8  subarrays  for  receive  and  a  continuous  35  by  35  element  transmit 
array. 
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Figure  17.  Configuration  of  MFAR  DSA 
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2.  Pattern  Synthesis  of  Shipboard  MFAR  DSA 

Based  on  the  synthesis  method  introduced  on  Section  A  of  this  chapter,  Figure  18 
shows  the  receive  pattern  from  coherent  signal  processing  of  16x8  uniform  subarrays. 
The  grating  lobes  exist  at  6  =  1 1.5°, 23°, 37°, 53°  (main  beam  not  scanned).  By 
multiplying  the  transmit  pattern  of  the  center  35  elements  as  shown  on  Figure  13,  the 
resulting  two-way  pattern  in  the  H  and  E-planes  is  given  in  Figures  19  and  20.  The  peak 
SLL  is  approximately  -45  dB.  The  grating  lobes  in  the  two-way  pattern  have  been 
reduced  to  a  maximum  of  —45  dB  below  the  two-way  main  beam  gain.  This  demonstrates 
the  effectiveness  of  the  2-way  pattern  approach. 


Figure  18.  Receive  H-plane  pattern  of  16  Hamming  weighted  subarrays  spaced  5 A. 
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Two-way  pattern  at*  =  90  deg  in  dB  Two-way  pattern  at*  =  0  deg  in  dB 
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Figure  19.  Two-way  H-plane  broadside  pattern  of  MFAR  DSA 


Figure  20.  Two-way  E-plane  broadside  pattern  of  MFAR  DSA 
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Figures  21  and  22  show  the  receive  patterns  when  the  beam  is  scanned  to 
Bs  =10°  and  60°  respectively.  The  nulling  effect  of  subarray  factor  on  the  configuration 
factor  is  unchanged  with  scan.  The  transmit  pattern  broadens  from  scanning  at  the  same 
rate  as  the  receive  pattern  broadens.  Figures  23  and  24  show  the  two-way  pattern 
scanning  to  0,  =10°  and  60°,  Figure  25  shows  the  two-way  pattern  of  multiple  beams 
scanned  at  increments  of  A0  =  5° . 


Figure  21.  Receive  pattern  of  16  subarrays  when  scanned  to  10°  from  H-plane 
broadside 
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Figure  22.  Receive  pattern  of  16  subarrays  scanned  to  60°  from  H-plane  broadside 


Figure  23.  Two-way  pattern  scanned  to  10°  from  broadside 
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Figure  24.  Two-way  pattern  scanned  to  60°  from  broadside 


Figure  25.  Multiple  beams  scanning  with  ^increments  5° 
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3.  Shipboard  HF  and  VHF  Distributed  Subarray  Antennas 

The  other  possible  application  of  DSA  onboard  a  ship  is  at  the  lower  frequencies. 
Lower  frequencies  have  better  performance  for  long-range  surveillance,  especially  on 
targets  with  small  RCS  because  HF  and  VHF  are  in  the  resonance  region  for  the  targets 
of  interest  (e.g.,  cruise  missiles).  However,  a  major  problem  is  how  to  increase  the 
resolution  from  the  limited  small  areas  on  a  ship.  By  processing  signals  coherently  from 
several  separate  subarrays,  the  half  power  beamwidth  can  approach  XI L,  where  L  can 
be  the  total  length  of  the  ship. 

Figure  26  is  the  side  view  of  an  AEGIS  cruiser,  with  some  possible  areas,  marked 
as  1  to  4,  on  which  subarrays  could  be  placed.  The  locations  are  detailed  in  Figure  27.  For 
simplicity,  it  is  assumed  that  these  four  subarrays  are  all  flat  and  rectangular  in  shape. 
The  number  of  elements  is  12,  6,  7  and  13,  respectively  (along  the  x-axis).  The  physical 
limitations  in  this  case  are  significant,  and  the  number  of  elements  is  so  few  that  a 
straight  forward  combination  of  the  four  areas  does  not  yield  any  reasonable  perfonnance 
at  VHF.  The  placement  of  elements  in  each  subarray  needs  to  be  optimized  by  the 
Minimax  method  in  Section  A  to  decrease  the  level  of  the  grating  lobes. 

By  fixing  the  position  of  the  edge  elements  in  each  subarray,  the  Minimax 
algorithm  finds  the  optimized  positions  for  minimum  peak  SLL  in  the  array  factor.  Since 
subarrays  1  and  2  and  subarrays  3  and  4  overlap  along  the  x-axis,  and  the  maximum 
possible  distance  is  between  subarrays  1  and  4,  there  are  only  two  possible  combinations 
of  subarrays  to  produce  a  low  SLL  pattern  for  the  requirement  of  high  resolution.  The 
first  combination  is  to  process  all  four  subarrays  coherently  as  a  monostatic  radar.  A  plot 
of  the  array  factor  and  the  equivalent  two-way  pattern  is  shown  in  the  Figures  28  and  29 
for  combination  1 . 

A  second  possible  combination  is  to  use  subarrays  1  through  3  as  the  receive 
pattern,  and  then  synthesize  a  transmit  pattern  for  subarray  4  by  the  equiripple  method 
discussed  in  the  Section  A.  A  plot  of  the  receive  pattern  (3  subarrays)  is  given  in  Figure 
30,  and  the  two-way  pattern  in  Figure  31.  The  optimized  element  positions  for  the  2 
combinations  are  shown  in  the  Tables  1  and  2. 


42 


Figure  26.  Side  view  of  the  Aegis  cruiser 
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Figure  27.  Possible  locations  of  VHF  subarrays  on  the  Aegis  cruiser 
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Figure  28.  Combination  1:  Pattern  of  the  4  subarrays  processed  coherently  to  form  a 
single  array  used  for  both  transmit  and  receive. 


Figure  29.  Combination  1:  Equivalent  2-way  pattern  of  the  monostatic  array. 
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Figure  30.  Combination  2:  Receive  pattern  using  3  subarrays 


Figure  31.  Combination  2:  Two-way  pattern  (3  receive  subarrays  and  1  transmit 
subarray). 
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Element  position 
in  A 

Subarray  1 

Subarray  2 

Subarray  3 

Subarray  4 

*i 

0 

5.1090 

10.7439 

11.5530 

*2 

0.0073 

6.0490 

10.7659 

11.5588 

x3 

0.5146 

6.5390 

11.2379 

12.0647 

*4 

2.0219 

7.0290 

11.7099 

12.5705 

*5 

2.5292 

7.5190 

12.1819 

13.9289 

*6 

3.0365 

7.6070 

12.6539 

14.5822 

x7 

3.5438 

13.5760 

14.6629 

x8 

4.0511 

14.9961 

x9 

4.5071 

15.0996 

*10 

5.0657 

16.6055 

Xu 

5.5730 

17.1113 

X12 

5.5800 

17.6171 

X» 

17.6230 

Table  1.  Optimized  locations  of  elements  in  the  4  subarrays  for  combination  1. 


Element  position 
in  A 

Subarray  1 

Subarray  2 

Subarray  3 

Subarray  4 

Xi 

0 

5.1090 

10.7439 

0 

x2 

0.0073 

6.0490 

10.7659 

0.5 

x3 

0.5146 

6.5390 

11.2379 

1 

*4 

2.0219 

7.0290 

11.7099 

1.5 

*5 

2.0464 

7.5190 

12.1819 

2 

*6 

2.8683 

7.6070 

13.5539 

2.5 

x7 

3.5438 

13.5760 

3 

x8 

3.6610 

3.5 

x9 

3.6706 

4 

*10 

4.4881 

4.5 

*11 

5.5730 

5 

*12 

5.5800 

5.5 

*13 

6 

Table  2.  Optimized  locations  of  elements  in  the  4  subarrays  in  combination  2. 
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Combination  1 ,  which  uses  all  4  subarrays,  gives  better  perfonnance  for  both  the 
peak  SLL  (-22.5  dB)  and  the  beamwidth  (2°)  for  the  round  trip  pattern  as  shown  on 
Figure  29.  However,  Combination  2  has  a  much  lower  average  SLL  than  Combination  1, 
and  the  individual  use  of  subarray  4  gives  more  flexibility  on  the  design  of  two-way 
pattern. 

These  two  distributed  VHF  subarray  designs  demonstrate  the  possibility  of 
maximizing  the  resolution  given  the  space  limitation  onboard  the  ship.  This  simple 
algorithm  of  pattern  design  does  not  show  the  full  potential  of  unequal  sized  or 
nonuniform  spacing  subarrays.  Other  suitable  algorithms,  like  the  Genetic  Algorithm 
(GA),  may  provide  better  solutions  for  these  types  of  problem  [Ref.  16]. 

4.  Calculation  of  Antenna  Parameters  on  MFAR  DSA  Design 

Since  the  DSA  has  a  large  spacing  between  subarrays,  the  simple  formulas  for 
antenna  parameters  need  modification  from  those  for  conventional  filled  arrays.  The 
following  paragraph  will  show  analytical  and  simulated  results  of  these  important 
parameters  of  the  MFAR  DSA.  The  pattern  parameters  will  be  examined  for  the  antenna 
design  described  in  Section  2. 

(a)  Two-way  beamwidth  between  first  nulls 

Since  the  result  of  the  two-way  pattern  is  the  product  of  transmit  and  receive 
patterns,  by  the  principle  of  pattern  multiplication,  the  two-way  BWFN  is  a  function  of 
both  patterns.  The  BWFN  of  the  receive  pattern,  from  Equation  (2),  is  mainly  dependent 
on  the  subarray  configuration  factor,  which  varies  as  sm{Ny/!2)  in  the  unifonn  DSA 
design.  The  BWFN  of  transmit  pattern,  however,  is  detennined  by  the  Chebyshev 
coefficients,  which  are  controlled  by  the  sidelobe  level,  number  of  elements,  interelement 
spacing  and  the  frequency.  The  transmit  BWFN  is  not  easy  to  calculate  analytically, 
especially  for  the  two-way  pattern.  Theoretically,  since  the  beamwidth  of  the  receive 
pattern  is  sharper  than  that  of  transmit  pattern,  the  shape  of  the  mainlobe  at  broadside  is 
determined  mainly  by  the  receive  pattern.  The  BWFN  of  the  receive  pattern  is 
approximately  3.2°  at  the  H-plane  broadside  measured  from  the  pattern  plot  on  Figure 
17. 
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(b)  Two-way  half  power  beamwidth 

Since  the  two-way  mainlobe  is  determined  mainly  by  the  receive  pattern,  the 
HPBW  of  the  two-way  pattern  can  be  approximately  detennined  by  the  receive  pattern 
also,  which  is 


\AFnorm(U)  f 


sin(^)  sin(^) 


M  sin(^)  N  sin(^) 


=  0.5. 


(14) 


It  is  assumed  unifonn  weighted  on  both  subarrays  and  elements.  Since  ME, ,  and 

the  phase  term  of  subarray  configuration  is  changing  much  faster  than  the  phase  of 
subarray  factor,  the  HPBW  is  detennined  mainly  by  the  subarray  configuration  factor. 


A  =  wavelength ;  D  =  aperture  wideth 


Type  of  distribution,  |z|  <  1 

Relative 

gain 

Half-power 

beamwidth, 

deg 

Intensity  of  first  sidelobe, 
dB  below  maximum  intensity 

Uniform;  A(z)  =  1 

1 

5 1  A/D 

13.2 

Cosine;  A(z)  —  cos"( 77^/2): 

n  =  0 

1 

51  AID 

13.2 

n  =  1 

0.810 

69A/D 

23 

n  =  2 

0.667 

S3A/D 

32 

n  =  3 

0.575 

95  A/D 

40 

n  =  4 

0.515 

1 1 1  A/D 

48 

Parabolic;  /t(z)  =  1  —  (1  -  Ajz2: 

A  =  1.0 

1 

51  A/D 

13.2 

A  =  0.8 

0.994 

53A/D 

15.8 

A  =  0.5 

0.970 

56A ID 

17.1 

A  =  0 

0.833 

bbXJD 

20.6 

Triangular;  A(z)  =  1  —  |;| 

0.75 

73A/D 

26.4 

Circular;  A(z)  =  V 1  -  z2 

0.865 

58.5  A/D 

17.6 

Cosine-squared  plus  pedestal; 

0.33  +  0.66  cos2  (irz/2) 

0.88 

63  A/D 

25.7 

0.08  +  0.92  cos2  (ttz/ 2),  Hamming 

0.74 

76.5A/D 

42.8 

Table  3.  Pattern  characteristics  produced  by  various  aperture  distributions  [from 
Ref.  10], 
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For  a  Hamming  window  distribution,  the  HPBW  is  approximately  76.5/1/  L  from 
Table  3.  For  /V  =  16,  D  =  5X,  L  =  15x5  +  2  =  11 X ,  by  the  parameters  of  MFAR  DSA 
design  in  Section  2,  the  broadside  HPBW  is  approximately  1°.  This  value  is  a  little  larger 
than  the  measured  value  of  0.9°  from  the  two-way  pattern  on  Figure  19  because  of  the 
multiplication  between  the  transmit  and  receive  pattern.  Compared  to  a  fully  populated 
uniform  array  (XI 2),  it  is  necessary  to  have  115  linear  uniform  weighted  elements  to 
have  the  same  resolution.  It  is  achieved  using  80  linear  elements  with  the  DSA  approach, 
(c)  Aperture  efficiency  and  directivity 

Since  the  receive  DSAs  are  Hamming  weighted  at  the  subarray  level,  the  taper 
efficiency  (receive  aperture  efficiency)  is  0.47  using  Equation  (3).  The  efficiency  of  the 
transmit  pattern  by  the  same  equation  is  0.6.  However  the  directivity  cannot  be  calculated 
by  the  equation  D  =  Am] A  cos#  /  X2 ,  because  grating  lobes  exist.  The  area  A  is  the  area 
occupied  by  array  elements. 

By  definition,  the  directivity  is  the  ratio  of  the  radiation  intensity  to  the  average 
radiation  intensity,  or  [Ref.  7] 

Air 

D(0, </>)  =  —,  (15) 

lia 


2k  n  „ 

sin  dddd(f)  is  the  beam  solid  angle,  and  E is  the  nonnalized 

0  0 

electric  field  intensity.  Neglecting  the  element  factor,  fs«0™  is  the  normalized  total  array 
factor.  By  numerical  integration  of  the  normalized  array  factor  in  spherical  coordinates, 
the  directivity  calculated  from  Equation  (15)  is  about  40  dB  for  a  uniform  weighted  DSA, 
which  is  almost  the  same  value  calculated  from  the  same  number  of  elements  ( 80  x  40 )  in 
contiguous  spacing  (XI 2).  Considering  the  tapering,  the  directivity  of  the  Hamming 
weighted  DSA  is  about  36.7  dB,  the  directivity  of  the  transmit  arrays  (35x35  elements) 
is  about  32.7  dB,  and  the  equivalent  one-way  antenna  gain  G  =  jG~Gr  -34.7dB.  This 

result  shows  that  the  gain  of  DSA  antenna  depends  directly  on  the  number  of  elements 
and  taper  efficiency. 


where  £1 


=  11 
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Table  4  summarizes  the  antenna  parameters  for  a  DSA  (receive  subarray)  design 
(N  =  16,  M  =  5,  Dx  =  5  A,  Dy  =  5A,  dx=0.5A,  dy=0.5A).  An  effective  area  of 

(0 . 5 A ) 2  is  assumed  for  each  element.  The  results  show  that  the  approximate  formula  for 
directivity  (D  =  47r7jAcos0J A1)  is  not  correct  in  the  DSA  application,  unless  the 
occupied  area  is  used. 


Approximate  formula 

(uniform  weight) 

Actual  computed 

(uniform  weight) 

Hamming  weight 

on  subarrays 

BWFN 

1.43° 

1.43° 

3.25° 

HPBW 

0.66° 

0.64° 

0.94° 

DIRECTIVITY 

45.5  dB  (total  area) 

40  dB  (active  area) 

40  dB 

36.7  dB 

MAX  SLL 

-13.2  dB 

-13.2  dB  (SLL) 

-3.8  dB  (grating) 

-39  dB  (SLL) 

-3.8  dB  (grating) 

Table  4.  Summary  of  antenna  pattern  parameters  for  MFAR  DSA. 


C.  SUMMARY 

This  chapter  has  presented  the  basic  theory  and  methods  of  subarraying.  There  are 
two  approaches  to  reducing  grating  lobes  due  to  the  large  spacing  between  subarrays  in  a 
periodic  DSA.  The  first  approach  is  placing  the  nulls  of  one  factor  ( AF  and  A Flx )  in  the 

direction  of  grating  lobes  of  the  other  factor  ( AFC ).  The  second  approach  is  to  use 

irregular  spaced  subarrays  to  reduce  the  peak  grating  lobes  by  redistributing  the  energy 
into  sidelobe  regions. 
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Two  possible  applications  of  DSA  onboard  a  ship  have  been  introduced  in  this 
chapter.  MFAR  DSA  can  be  used  at  higher  frequencies  (X  or  C  band)  to  reduce  the  cost 
and  weight  of  modern  radar.  At  lower  frequencies  (HF  or  VHF  band),  the  application  of 
DSA  can  increase  the  resolution  obtained  from  utilizing  the  limited  small  areas  on  a  ship. 
This  is  important  for  detecting  and  tracking  targets  with  small  RCS  like  cruise  missiles. 

The  calculation  of  antenna  parameters  of  the  DSA  design  shows  the  tradeoff  in 
performance  on  the  directivity  and  beamwidth.  The  directivity  depends  on  the  number  of 
elements  and  tapering  only;  the  beamwidth  decreases  as  the  subarray  spacing  increases. 
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IV.  CONCLUSION 


A.  SUMMARY 

The  concept  of  distributed  subarray  antennas  has  been  proposed  for  both  the 
MFAR  and  VHF  applications.  By  combining  distributed  subarrays  on  the  available  areas 
of  a  constrained  platform,  the  MFAR  or  VHF  DSA  can  achieve  the  maximum  resolution 
(aside  from  synthetic  aperture  approaches)  and  potential  reductions  in  cost  and 
complexity.  The  two-way  pattern  design  of  a  DSA  effectively  suppressed  the  undesired 
grating  lobes  by  the  design  of  separate  transmit  and  receive  patterns.  From  the  pattern 
multiplication  principle,  the  grating  lobes  in  the  subarray  configuration  pattern  ( AFc ) 

have  been  suppressed  by  the  design  of  subarray  pattern  ( AFs )  and  transmit  pattern 

( afTx ). 

The  design  examples  of  shipboard  MFAR  have  shown  that  the  HPBW  decrease 
(increased  resolution)  can  be  achieved  by  spreading  a  fewer  number  of  elements  over  a 
longer  baseline  and  then  suppressing  the  grating  lobes  with  other  pattern  factors.  The 
nulling  effect  of  the  subarray  factor  does  not  change  as  the  beam  scans  away  from 
broadside.  The  limitations  and  advantages  of  DSA  have  been  mentioned  previously 
throughout  the  thesis,  and  are  summarized  here. 

1.  Advantages  of  the  MFAR  DSA 
a.  Support  of  Wideband  ADBF 

The  collection  of  elements  is  steered  in  phase  at  the  element  level,  and  a 
beamforming  network  combines  the  element  into  subarrays.  The  subarrays  are  then 
steered  via  photonic  time  delay  devices.  This  provides  the  necessary  wideband  beam 
steering  and  ADBF  capabilities  at  a  reasonable  cost. 
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b.  Support  of  Multi-functionality 

Since  each  subarray  is  physically  separated,  a  separate  function  can  be 
assigned  each  subarray  more  directly.  The  design  of  subarrays  can  be  more  adaptive  to 
the  main  function  from  the  design  stage,  which  decreases  the  effects  of  compromises  on 
performance. 

c.  Low  Complexity  of  Manufacturing  and  Computing  with  Less  Control 
Elements 

Although  some  compromises  between  performance  and  periodic  spacing 
have  been  made,  regularity  makes  the  implementation  of  subarrays  more  realistic  and 
practical  for  the  requirements  of  modern  radar.  However,  irregular  placement  of  elements 
is  not  a  huge  disadvantage.  Randomly  thinned  arrays  have  been  used  for  radars  in  the  past 
[Ref.  16]. 

d.Achieving  High  Resolution  with  Less  Space  Limitation 
The  required  angular  accuracy  and  range  resolution  can  be  achieved  with 
little  limitation  on  platfonn  space.  Little  perturbations  in  subarray  locations  do  not  have 
much  effect  on  the  array  pattern.  The  arrangement  and  dimensions  of  the  subarrays  can 
be  adjusted  to  the  shape  of  platfonn. 

e.  Possibility  of  Multi  Band  Shared  Apertures 

Since  the  spacing  between  the  subarrays  is  large,  it  is  possible  to  insert 
other  low  frequency  elements  to  share  the  same  aperture  space.  For  example,  if  the 
subarrays  structure  is  designed  for  X-band  frequencies,  then  an  L-band  array  element  can 
be  inserted  between  subarrays  with  half  wavelength  spacing.  Figure  32  shows  a  possible 
scheme  for  a  X  and  L  band  array. 
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X  band  subarray 


L  band 
element 


L  band 
element 


X  band  subarray 


PHASE 

SHIFTERS 

SUBARRAYS 

IN/OUT 


LOW  NOISE 
AMPLIFIER 


TRUE 

TIME 

DELAY 


ARRAY  IN/OUT 


L/X  BAND 


Figure  32.  Possible  beam  control  scheme  of  a  dual-band  DSA  [after  Ref  17] 


2.  Limitations  of  the  MFAR  DSA  Design 

a.  Less  Efficiency  on  Surveillance  Due  to  the  Narrow  Beamwidth 
of  Transmit  Pattern 

If  the  beamwidth  of  the  transmit  pattern  is  tied  to  the  synthesis  of  the  two- 
way  pattern,  the  resulting  beam  may  not  provide  a  large  enough  field  of  view  for 
searching.  This  disadvantage,  however,  can  be  corrected  by  dividing  the  transmit  pattern 
into  searching  and  tracking  modes.  In  the  search  mode,  the  center  transmit/receive 
module  operates  alone  like  any  other  element-based  phased  array.  In  the  tracking  mode, 
the  center  transmit  module  is  a  high  contrast  pencil  beam,  and  the  resultant  two-way 
pattern  is  the  multiplication  of  transmit  and  receive  pattern  synthesized  for  the  DSA.  A 
dual-band  shared  aperture  is  another  possible  approach  for  separating  the  track  and 
surveillance  into  different  combinations;  for  examples,  long-range  surveillance  (L  band), 
mid-range  surveillance  (X  band)  and  tracking  (X  band  DSA). 
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b.  Lower  Beam  Efficiency  Due  to  the  Grating  Lobes 

From  the  calculations  in  Chapter  III,  the  antenna  directivity  depends  only  on 
the  area  occupied  by  elements.  There  is  a  reduction  in  beam  efficiency  resulting  from  the 
grating  lobes,  and  any  sidelobe  taper  will  lower  the  efficiency  further  (typically  about  3 
dB).  However,  this  can  be  improved  at  the  expense  of  more  subarrays  or  closer  spacing 
between  them 

c.  Fewer  Degrees  of  Freedom  on  Pattern  Synthesis 

Since  the  number  of  individual  control  elements  is  decreased  by  the  subarray 
design,  the  degrees  of  freedom  (DOF)  in  terms  of  ADBF  capability  is  decreased 
substantially.  The  limitations  resulting  from  the  synthesis  of  transmit  and  receive  patterns 
also  restrict  the  implementation  of  DSA  in  the  environment  of  high  clutter  because  of 
lower  sensitivity  and  less  capability  for  nulling  interferences. 

d.  More  Complexity  of  Task  Scheduling 

Having  a  multifunction  radar  complicates  the  task  scheduling  of  a  single 
antenna.  The  two-way  pattern  approach  introduces  more  constraints  on  the  scheduling 
tasks.  The  time  budget  i.e.,  the  allocation  of  radar  time  to  different  tasks  is  dependent  on 
the  radar  antenna  parameters.  Since  the  main  tasks  of  MFAR  are  search,  track  and  other 
auxiliary  functions,  the  strategy  of  how  to  execute  of  all  tasks  in  the  best  possible  way  is 
very  complicated. 


B.  POSSIBLE  FURTHER  RESEARCH  TOPICS 

There  are  two  primary  directions  of  research  on  suppressing  the  grating  lobes  in 
the  pattern  of  widely  spaced  subarrays,  as  discussed  below. 

1.  Filtering  Approach 

There  are  many  filtering  techniques  already  developed  in  the  area  of  digital  signal 
processing.  If  the  main  objective  is  to  suppress  the  grating  lobes  from  self-interference, 
the  adaptive  notch  filter  (ANF)  or  space-time  domain  adaptive  processing  (STAP)  might 
be  the  possible  research  area  for  the  implementation  of  DSA. 
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Adaptive  notch  filters  are  widely  used  in  many  signal-processing  applications  to 
extract  and  trace  the  narrow-band  noise.  The  basic  principle  is  to  produce  the  adaptable 
frequency  response  in  both  the  time  and  space  domain:  zero  at  the  specified  spatial  angle, 
one  otherwise.  Figure  33  shows  the  example  of  frequency  response  on  nulling  the  angles 
±15°  and  ±37°,  everywhere  else  unchanged.  The  depth  and  width  of  nulls  should  be 
adjusted  according  the  grating  conditions.  But  this  must  be  implemented  in  both  the  space 
and  time  domain,  which  introduces  the  applications  of  the  STAP  algorithm  [Ref.  18].  For 
a  DSA  the  notches  would  be  placed  at  grating  lobe  locations. 

STAP  processes  signals  in  the  spatial  and  time  domains.  The  time  domain 
includes  both  slow-time  (pulse  repetition  interval)  and  fast-time  (range  cell).  Therefore, 
the  utilization  of  phase  and  amplitude  weighting  as  a  means  to  achieve  a  desired  steering 
or  nulling  direction  must  consider  the  spatial,  slow-time  and  fast-time  factors. 

2.  Digital  Arrays 

Unlike  an  analog  beamforming  network,  digital  beamfonning  arrays  digitize 
received  signals  at  the  element  level,  then  process  these  signals  in  a  digital  processor  to 
form  the  desired  beam  and  frequency  response  (in  this  case,  nulling  in  the  grating 
direction).  By  doing  so,  the  total  information  available  at  the  aperture  has  been  properly 
preserved  and  can  be  manipulated  indefinitely  without  introducing  further  error  (other 
than  computer  roundoff  error). 

Adaptive  digital  beamforming  at  the  element  level  can  reject  interference  and  at 
the  same  time  steer  a  main  lobe  in  the  direction  of  a  desired  signal.  Therefore,  any 
number  of  beams  can  be  fonned  or  the  signal  can  be  rejected  according  to  a  selected 
algorithm. 
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Figure  33.  Adaptive  nulling  at  the  angle  ±15°,  ±37°  by  ANF. 
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APPENDIX  A:  GLOSSARY  OF  TERMINOLOGY 


MFAR 

ADBF 

LPI 

MFSAR 

dB 

mrad 

MHz 

GHz 

MESAR 

TWS 

AF 

AFS 

AFC 

afTx 

afKx 

EF 

BWFN 

HPBW 

X 

SNR 

RCS 


Multi-function  Array  Radar 
Adaptive  Digital  Beamforming 
Low  Probability  of  Intercept 
Multifunction  Subarray  Radar 
Decibels 

milli  radian  ( 1 0  radian) 

Megahertz  (106  cycles/second) 

Gigahertz  ( 109  cycles/second) 

Multifunction  Electronically  Scanned  Adaptive  Radar 
Track- While-Scan 
Array  Factor 
Subarray  Factor 

Subarray  Configuration  Factor 
Transmit  Array  Factor 
Receive  Array  Factor 
Element  Factor 

Beamwidth  Between  First  Nulls 
Half  Power  Beamwidth 
Lambda  —  wavelength 
Signal-to-Noise  Ratio 
Radar  Cross  Section 
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SLL 

Sidelobe  Level 

DSA 

Distributed  Subarray  Antenna 

FIR 

Finite  Impulse  Response 

WLS 

Weighted  Least  Squares  algorithm 

Tx 

Transmitter 

Rx 

Receiver 

HF 

High  Frequency  (3  ~  30  MHz) 

VHF 

Very  High  Frequency  (30  ~  300  MHz) 

UHF 

Ultra  High  Frequency  (300  ~  3000  MHz) 

DOF 

Degree  Of  Freedom 

ANF 

Adaptive  Notch  Filter 

STAP 

Space  Time  Adaptive  Processing 
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APPENDIX  B:  MATLAB  CODE  LISTING 


The  major  MATLAB  Code  listing  for  the  pattern  plot  as  the  following. 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  DSA.m 

%  Pattern  calculation  and  plot  for  Figure  8,13,14,18-25 
%  Quasi-monostatic  subarrays  AF 

%  XMTR  AF  has  35x35  Remez-weighted  elements  spacing  0.5  wavelength 
%  RCVR  is  16x8  subarrays  with  5x5  elements  in  each  subarray 
%  Each  subarray  spacing  5  wavelengths  in  x&y  axis  with  Flamming  weights 
%  Display  the  aperture  efficiency  also 
%  resolution  of  plot  is  quarter  degree 

clear  all; 

llx=0:5:75;  %  16  Subarray  spaced  (Lambda)  in  x-axis 
lly=0:5:35;  %  8  Subarray  spaced  in  y-axis 

dx=0.5;  Mx=5;  %  element  spacing  and  number  in  subarray  along  x-axis 
dy — 0.5 ;  My=5;  %  element  spacing  and  number  in  subarray  along  y-axis 

%  Scan  angle  input 

%  Figure  21-25  need  to  change  angle  accordingly 

thets=0; 

phis=0; 

us=sin(thets*pi/l  80)*cos(phis*pi/l  80); 
vs=sin(thets*pi/180)*sin(phis*pi/180); 

thet=linspace(-90,90,720);  %  x,y  direction  cosine 
theta=thet.'; 

phi=linspace(0, 180,720); 
u=sin(theta*pi/l  80)*cos(phi*pi/l  80); 
v=sin(theta*pi/180)*sin(phi*pi/180); 

U=u-us; 

V=v-vs; 

%  Subarray  configuration  AF 
%  Weighting  Function  by  Flamming  window 
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wx=zeros(l,length(llx));  %  weights  in  x-axis 
for  n=l  :length(llx) 

xn=(n- 1  )/(length(llx)- 1  )*2- 1 ; 
wx(n)=0.92*abs(cos(xn*pi/2))A2+0.08; 
end 

wy=zeros(l,length(lly));  %  weights  in  y-axis 
for  n=l  :length(lly) 

yn=(n-  l)/(length(lly)- 1  )*2- 1 ; 
wy(n)=0.92*abs(cos(yn*pi/2))A2+0.08; 
end 

C  AF  x=zeros  ( 72 0 ) ; 
for  lx=l:length(llx) 
psilx=2*pi*llx(lx)*U; 
CAFx=CAFx+wx(lx)*exp(J*psilx); 
end 

CAFy=zeros(720); 
for  ly=l  dength(lly) 
psily=2*pi*lly(ly)*V; 
CAFy=CAFy+wy(ly)*exp(j*psily); 
end 

CAF=CAFx.*CAFy/sum(wx)/sum(wy); 

%  Single  subarray  AF 
SAFx=zeros(720); 
for  sx=  1  :Mx 

psix=2*pi*dx*U; 

SAFx=SAFx+exp(j*(sx-l)*psix); 

end 

SAFy=zeros(720); 
for  sy=l:My 
psiy=2*pi*dy*V; 
SAFy=SAFy+exp(j*(sy-l)*psiy); 
end 

SAF=SAFx.  *  S  AFy/Mx/My ; 
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%  XTMR  array  design  to  pass  mainlobe,  others  with  45  db  attenuation 
%  First  calculation  the  equal  ripple  coefficients 

%  use  for  obtaining  9  deg  mainlobes,  every  other  angle  response  -40  db 

px=35;  %  Element  numbers  in  x-axis 

fs=90;  %  Array  factor  half  space 

rp=l ;  %  Passband  ripple  in  db 

rs=-50;  %  Stopband  desired  attenuation  in  db 

fc=[0  4.5];  %  Cutoff  spatial  freq 

a=[l  0];  %  Desired  amplitude 

dev=[(10A(rp/20)-l)/(10A(rp/20)+l)  10A(rs/20)];  %  dB  to  numerical 
[nx,fOx,aOx,Wx]=remezord(fc,a,dev,fs); 
bx=remez(px-l,fOx,aOx,Wx);  %  Using  fixed  number  elements 
bx=abs(bx(l:length(bx)));  %  Elements  weighting 
Bx=bx/max(bx);  %  Normalized  weights 

py=35;  %  element  #  in  y-axis 

fs=90;  %  Array  factor  half  space 

rp=l ;  %  Passband  ripple  in  db 

rs=-50;  %  Stopband  attenuation  in  db 

fC=[0  4.5];  %  Cutoff  spatial  freq(angle) 

a=[l  0];  %  Desired  amplitude 

dev=[(  1 0A(rp/20)- 1 )/( 1 0A(rp/20)+ 1 )  1 0A(rs/20)] ; 

[ny,fOy,aOy,Wy]=remezord(fC,a,dev,fs); 

by=remez(py-l,fOy,aOy,Wy);  %  Using  fixed  number  elements 

by=abs(by(l:length(by)));  %  elements  weighting 

By=by/max(by);  %  Normalized  weights 

%  Array  factor  in  square  XTMR 
XAFx=zeros(720); 
for  Xx=l  dength(Bx) 
psiXx=2*pi*dx*U; 

XAFx=XAFx+Bx(Xx)*exp(j*(Xx-l)*psiXx); 

end 

XAFy=zeros(720); 
for  Xy=l  dength(By) 
psiXy=2*pi*dy*  V ; 
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XAFy=XAF  y+By  (Xy )  *  exp(j  *  (Xy- 1 )  *psiXy ) ; 
end 

%  calculate  RCVR  aperture  efficiency 
Rsl=0;  Rs2=0; 
for  il=l:length(llx) 
for  i2=l:length(lly) 

Rsl=Rsl+abs(wx(il)*wy(i2));  Rs2=Rs2+abs(wx(il)*wy(i2))A2; 
end 
end 

eta=RslA2/length(llx)/length(lly)/Rs2; 

%  calculate  XTMR  aperture  efficiency 
Xsl=0;  Xs2=0; 
for  il=l:length(Bx) 
for  i2=l  :length(By) 

Xsl=Xsl+abs(Bx(il)*By(i2));  Xs2=Xs2+abs(Bx(il)*By(i2))A2; 
end 

end 

ETA=Xs  1  A2/length(Bx)/length(By)/Xs2; 
disp(['RCVR  aperture  efficiency:  ',num2str(eta)]) 
disp(['XTMR  aperture  efficiency:  ',num2str(ETA)]) 

%  XTMR  AF 

XAF=(XAFx.*XAFy)/sum(Bx)/sum(By); 

dbXAF=20*logl0(abs(XAF)); 

%  RCVR  AF 

RAF=SAF.*CAF; 

dbRAF=20*logl0(abs(RAF)); 

%  Two-way  pattern 
BAF=XAF.*RAF; 

dbBAF=20*logl0(abs(BAF));  %  2-way  pattern  in  dB 
dbSAF=20*logl0(abs(SAF));  %  Subarray  pattern  in  dB 
dbCAF=20*logl0(abs(CAF));  %  Subarray  configuration  pattern  in  dB 

figure(l)  %  Figure  8 

subplot(3,l,l),  plot(theta(  1:720),  dbRAF((l:720),phis+l)),axis([-90  90  -60  0]), 
grid  on,  title('  Receive  pattern  of  16  subarrays  spacing  5\lambda') 
subplot(3,l,2),  plot(theta(  1:720),  dbSAF((l:720),phis+l)),axis([-90  90  -60  0]), 
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grid  on,  title('  Pattern  of  single  subarray  which  is  composed  by  5  elements') 
ylabel('Relative  pattern  at  \phi=0\circ  in  dB') 

subplot(3,l,3),  plot(theta(  1:720),  dbCAF((l:720),phis+l)),axis([-90  90  -60  0]), 
grid  on,  title)'  Pattern  of  16  subarray  configuration') 
xlabel('\theta  in  degree') 

figure(2) 

subplot(3,l,l),  plot(theta)  1:720),  dbBAF((l:720),phis+l)),axis([-90  90  -100  0]), 
grid  on,  title)'  Two  Way  pattern  of  16  subarrays  spacing  5\lambda') 
subplot(3,l,2),  plot(theta)  1:720),  dbRAF((l:720),phis+l)),axis([-90  90  -100  0]), 
grid  on,  title)'  Receive  pattern  of  16  subarrays  spacing  5\lambda') 
ylabel('Relative  pattern  at  \phi=0\circ  in  dB') 

subplot(3,l,3),  plot(theta)  1:720),  dbXAF((l:720),phis+l)),axis([-90  90  -100  0]), 
grid  on,  title)'  Transmit  pattern  of  35  elements') 
xlabel('\theta  in  degree') 

figure(3)  %  Figure  18 

%  Receive  pattern  of  16  subarrays  spaced  5  wavelength 
plot) theta)  1:720),  dbRAF(l:720,phis+l)),axis([-90  90  -60  0]), 
grid  on, 

ylabel('Receive  pattern  at  \phi  =  0  deg  in  dB') 
xlabel('\theta  in  degrees') 

figure(4)  %  Figure  13 
%  Transmit  pattern  of  35  elements 

plot) theta)  1:720),  dbXAF(l:720,phis+l)),axis([-90  90  -80  0]), 
grid  on, 

xlabel('\theta  in  degrees') 

ylabel('Transmit  pattern  at  \phi  =  0  deg  in  dB') 

figure(5)  %  Figure  19 

%  plot  Two-way  pattern  of  16x8  subarrays  spaced  5  wavelength  at  Fl-plane 
plot) theta)  1:720),  dbBAF(l:720,l)),axis([-90  90  -80  0]), 
grid  on,ylabel('Two-way  pattern  at  \phi  =  0  deg  in  dB'), 
xlabel('\theta  in  degrees') 

figure(6)  %  Figure  20 
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%  plot  Two-way  pattern  of  16x8  subarrays  spaced  5  wavelength  at  E-plane 
plot(theta(  1:720),  dbBAF(l:720,361)),axis([-90  90  -80  0]), 
grid  on,ylabel('Two-way  pattern  at  \phi  =  90  deg  in  dB'), 
xlabel('\theta  in  degrees') 

%  End  of  DSA.m 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  IrregSubarray.m 

%  AF  of  VHF  DSA  with  freq  100  MHz 

%  Looking  for  the  locations  of  38  optimized  element  in  4  subarrays 
%  Find  element  optimized  location  and  first  null  by  myfunlO.m 
%  All  position  in  wavelengths 
%  Used  for  Figure  28,  29  pattern  plot 

clear  all; 

u0=0.05;  %  Initial  guess  of  beamwidth 
ul=l; 

%  element  locations  of  first  subarray 

s0=0;  %  Lower  boundary  of  elements  in  1st  subarray 

si  1=5.58;  %  Upper  boundary  of  elements 

for  n=l:10 

sl(n)=0.5073*n; 

end 

%  Initial  guess  and  boundary  of  12  element  locations 
sl0=cat(2,s0,sl,sl  1); 
sll=cat(2,s0,sl-0.5,sl  1); 
slu=cat(2,s0,sl+0.5,sl  1); 
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%  element  locations  of  2nd  subarray 
s21=5.109;  %  First  element  position 
s26=7.607;  %  Last  element  position 
for  n=l:4 

s2(n)=s21+n*0.49; 

end 

%  Initial  guess  and  boundary  of  6  element  locations 
s20=cat(2,s2 1  ,s2,s26); 
s21=cat(2,s21,s2-0.45,s26); 
s2u=cat(2,s2 1  ,s2+0.45,s26); 

%  element  locations  of  3rd  subarray 
s3 1=1 0.7439;  %  First  element  position 
s37=13.576;  %  Last  element  position 
for  n=l:5 

s3(n)=s31+n*0.472; 

end 

%  Initial  guess  and  boundary  of  7  element  locations 
s30=cat(2,s31,s3,s37); 
s31=cat(2,s31,s3-0.45,s37); 
s3u=cat(2,s3 1  ,s3+0.45,s37); 

%  element  locations  of  4th  subarray 
s41=l  1.553;  %  First  element  position 
s4 1 3=  1 7.623 ;  %  Last  element  position 
forn=l:ll 

s4(n)=s41+n*0. 50583; 
end 

%  Initial  guess  and  boundary  of  13  element  locations 
s40=cat(2,s4 1,s4,s4 1 3); 
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s41=cat(2,s41,s4-0.5,s413); 

s4u=cat(2,s41,s4+0.5,s413); 

%  Combination  of  all  element  positions 

Ilb=cat(2,u0,sll,s21,s31,s41);  %  lower  bound  of  element  and  the  first  null  position 
Iub=cat(2,ul,slu,s2u,s3u,s4u);  %  upper  bound  of  element  and  the  first  null  position 
10=cat(2,u0,sl0,s20,s30,s40);  %  Intial  guess  of  element  and  the  first  null 

options=optimset('MinAbsMax',  1 ); 

[l,fVal,  maxfval,  exitflag,  output]=fminimax(@myfunlO, 10, [],[],[],[], lib, lub,[], options) 

%  Result  pattern  calculation 
%  Scan  angle  input  in  degrees 
thets=0;  phis=0; 

us=sin(thets*pi/180)*cos(phis*pi/180); 

vs=sin(thets*pi/180)*sin(phis*pi/180); 

%  x,y  direction  cosine 
thet=linspace(-90, 89. 75,720); 
theta=thet.'; 

phi=linspace(0, 179.75,720); 
u=sin(theta*pi/180)*cos(phi*pi/180); 
v=sin(theta*pi/l  80)*sin(phi*pi/l  80); 

U=u-us;  V=v-vs; 

%  AF  of  optimized  location 
L  AF  x=zero  s  (72 0 ) ; 
for  lx=l:length(ll) 
psilx=2*pi*ll(lx)*U; 

LAFx=LAFx+exp(j*psilx); 

end 

dbLAF=20*logl0(abs(LAFx)/length(ll)); 
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db2LF=40*logl0(abs(LAFx)/length(ll));  %  Equivalent  2way  pattern 


%  Pattern  plot  of  optimized  position  used  on  Fig.  28 
figure(l) 

plot( theta,  dbLAF(:,4*phis+l)),axis([-90  90  -60  0]), 
grid  on, 

xlabel('\theta  in  degrees') 
ylabel('Relative  pattern  in  dB') 

%  Equvalent  2-way  pattern  used  on  Fig.  29 
figure(2) 

plot( theta,  db2LF(:,4*phis+l)),axis([-90  90  -60  0]), 
grid  on, 

xlabel('\theta  in  degrees') 
ylabel('Relative  pattern  in  dB') 

%  end  of  IrregSubarray.m 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  myfunlO.m 

%  Function  for  optimization  location  of  subarrays 
%  Function  called  by  the  IrregSubarray.m 

%  Minimaxing  the  peak  sidelobe  level  in  the  visible  region[-pi/2,pi/2] 

function  F=myfunl0(l) 

%  Array  factor  of  same  subarray  with  irregular  spacing 
%  Computation  for  U  axis  only 
u0=linspace(l(  l),0.5*pi,5000); 


alpha=2*l; 
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F 1  =zeros(  1  ,length(uO)); 
for  n=2:39 

Fl=Fl+exp(j*alpha(n)*uO); 

end 

F=abs(Fl); 

%  end  of  myfunlO.m 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
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